Spectral diffusion on ultralong time scales in low temperature glasses 
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A dynamical theory is constructed to describe spectral diffusion in glasses in the temperature 
range near 1 Kelvin on long time scales. The theory invokes interacting tunneling centers (TLS) 
which provide an excess contribution to the spectral hole width which qualitatively accounts for the 
deviation from standard logarithmic line broadening observed by Maier et al. [Phys. Rev. Lett. 
76, 2085 (1996)]. Alternative explanation schemes of the nonlogarithmic line broadening, avoiding 
interacting TLS, are discussed. We devise experimental tests which could be used to access the 
validity of the proposed theories. 



PACS numbers: 61.43.FS, 78.40.-q 
I. INTRODUCTION 

At low temperatures the thermal, acoustic and opti- 
cal behavior of glasses differs significantly from that of 
crystalline solids. As examples note the linear temper- 
ature dependence of the specific heat, the pronounced 
absorption of sound even below 1 K, and the anomalous 
broadening of spectral holes in the homogeneous line of 
chromophore molecules. R is known now for more than 
twenty years that phonons cannot account for these ob- 
servations; instead localized low-energy excitations in the 
glass are needed. Since the introduction of the stan- 
dard tunneling model (STM) by Anderson, Halperin and 
VarmaEJ and PhillipsB, and the experimental observation 
of saturabilityjXif ultrasound by Hunklinger et alfl and 
Golding et al.Q, it is widely accepted that these low- 
energy excitations are two-level tunneling systems (TLS). 
Denoting the left and right ground state of the double- 
well potential by \L) and respectively, the Haniilto- 
nian reads in a Pauli-spin representation 



H 



(1) 



where A/?i is the tunneling frequency, e the asymmetry 
energy and u^, = \L){L\ - \R){R\. In the STM the inter- 
action between TLS is neglected, and it is assumed that 
the tunneling parameters A, e are random variables with 
distribution 



P(i)(e, A)dedA 



dedA, A > A„ 



(2) 



(with Po ~ 0.6 X 10-*^ ToT^ in PMMA, a polymer 
glass). The ensuing constant distribution for the TLS- 
energy splitting 



VA2 + e2 



(3) 



explains in particular the linear specific heat. Including 
relaxation of the TLS via the one-phonon process with 
rate 



R = 



R-n 



where [x = E/2kBT) 

R^^^{E) = aT^x^ coihx 



and 



(4) 



(5) 



(6) 



the model describes most acoustic and optical experi- 
ments in glasses satisfactorily. Here, 7 is the deformation 
potential energy of the TLS-phonon coupling, g the mass 
density of the glass, and v the sound velocity. 

Recent experiments report a systematic disagreement 
with the STM. An example is the attenuation of sound 
below 100 mK. The STM predicts a T'^-increaae in con- 
tradiction with the experimental r*^^~^^-law.Eli3 Devia- 
tions from the predicted STM behavior was also ob- 
served by Maier, Kharlamov, and|Jlaa,rer in low tem- 
perature hole-burning experiments Jl^'tHI They performed 
hole-burning of a chromophore embedded in PMMA at 
temperatures around 1 K up to extremely long times 
(from 10 s to 10 days). They found a logarithmic time- 
dependence with a crossover to an algebraic behavior af- 
ter about 3 h. Though the log t-behavior is in agreement 
with the STM, the algebraic behavior is not. The authors 
could fit their data with an ad hoc ansatz 



P{e,A)=Po 



1 

A 



A 
A2 



A ~ constant. 



(7) 



for the TLS-parameter distribution function and the as- 
sumption that relaxation occurs via the one-phonon pro- 
cess with the rate (^-(H). They motivate the distribution 
function (^ by recent publications focusing on the inteiu 
action of TLS in glasses. EjEj Indeed, Burin and KaganL^ 
have shown, following earlier ideas of Yu and Leggettti 
that pairs of interacting TLS do provide a means of con- 
stituting a distribution like the second term in Eq. (^ for 
certain excitations in their energy spectrum. They called 



TLS which are distributed according to the first term in 
(0) primary TLS, and those distributed according to the 
second term in i.e., pairs of primary TLS, secondary 
TLS. The nice feature of their theory is that the distribu- 
tion of secondary TLS [second term in Eq. (j^)] is derived 
from primary TLS which are distributed according to the 
STM [first term in Eq. (0)]. In that sense, Burin and Ka- 
gan's theory stays within the framework of the STM, and 
includes only excitations which have not been considered 
in the traditional treatment. So far these ideas have been 
worked out very qualitatively, and with emphasis on ex- 
periments in the millikelvin regime. Hence, the question 
arises whether they apply to the experiment of Maier 
et al., i.e., for relaxation processes at 1 K on the time 
scale between hours and days, or whether an extension 
of the STM — which would be as phenomenological as the 
STM — has to be found in order to understand this exper- 
iment. It is the purpose of this paper to address this issue 
by presenting a detailed model that includes TLS-TLS 
coupling in the manner proposed by Burin and Kagan. 

The paper is laid out as foll ows : in Sec. II, we analyze 
the experiment of Maier et al.tlHia and show the achieve- 
ments and failures of Burin and Kagan's approach for 
experiments in the Kelvin regime; in Sec. Ill and IV, we 
propose a specific model that combines interacting TLS 
with strong-coupling effects between TLS and phoncms 
in the framework of the theory of Kassner and SilbeyllZl; 
in Sec. V, we compare the predictions of the new model 
with the hole-burning data in Ref. and also discuss 
alternative explanation schemes, which comprise an ex- 
tension of the STM, and compare with equivalent hole- 
burning measurements in proteins; finally, in Sec. VI, we 
discuss our results and conclude with a short summary 
in Sec. VII. The mathematical details are relegated to 
three Appendices in order not to obscure the basic ideas. 



II. HOLE-BURNING AT ULTRALONG TIMES 
AND THE BURIN-KAGAN THEORY 



according to the following equation: 



(9) 



2 fco P f 

r(t).-(C)y^ dEsec^'^l dr-PiE,r) 



l-e 



-rR„ 



(10) 



Here, (C) is the chromophore-TLS coupling strength. 
The distribution P{E, r) directly follows from Eq. (Q) 
by using (|)-(|) and (|). With this, we find 



m = T{to) + — (C) Po 



ksT log^ + AV^(Vt-yj% 
to V 



(11) 



with 



a = 9 X 10^ K-^s"^ 



(12) 

in PMMA.0 To compare with the experiment, we have 
included in r{to) processes which are faster than the 
shortest experimental time to, for which a hole broad- 
ening can be determined. In Ref. 12(a) a fit to the 
experimental data provided 



IQ-'^ K 



(13) 



at 1 and 0.5 K with an error of approximately 10%. Fur- 
thermore, the experimental observation of spectral diffu- 
sion up to imax ~ 10^ S SUggCStS 



IQ-^ K 



(14) 



< 



1/Rn 



according to the relation 
(2A:BT/A„,in)2(l/ar3). 

Very recently. Burin and Kaganll3 showed that weak 
TLS-TLS interaction provides ultralow energy excita- 
tions. They added to the Hamiltonian (^ an interaction 
term 



i^TLS-TLS = -T X! ■^iJ^^^^^ 



In Ref. |11],|1^ photo-chemical hole- burning in PMMA at 

1 and 0.5 K has been performed for extremely long times 

imax = 10^ s. The authors found a log i-dependence with 

a crossover to an algebraic behavior after approximately 

3 h. The crossover shifts 1 order of magnitude in time , , ■ 

r 1 n4 J- 1 n3 --c j-u j- i ■ ■ j i, where the interaction energy 

from 10 to 10 s, if the temperature is increased by a 

factor of 2. 

A theoretical description of spectraLdiffusion in glasses 
was provided by EIu and WalkerQ and Black, and 
HalperinEl. ReineckeB, and later Bai and FayerE2l, ex- 
tended their results to optical experiments. Based on 
this work the dependence of the hole width V{t) on the 
waiting time t is determined by 



J i-'i 



(15) 



(16) 



PiE,r) 



Pq 



2vr 



1 



A 



(8) 



where 



falls off with distance |rj — r^l in a manner typical for a 
dipole-dipole interaction. They assumed that the angular 
average of the coupling is zero, 

= 0, (17) 

and that the TLS-TLS coupling is weak, i.e., 

PoUo « 1, (18) 



2 



where Uq is set by the variance 

(4> = ul 



(19) 



Such an interaction could be mediated by virtual phonon 
exchange (elastic coupling) or virtual photon exchange 
(electrostatic coupling) between the TLS. In the former 
case the energy scale of interaction (|9|) is easily found 
to be 



Uo 



7 



(20) 



which, iadeed, provides PqUq ^ 1 for all glasses — in 
PMMAta: [/q w 1 eV and PqUq w 10"'*. Based 
on the smallness of PqUq, Burin and Kagan proposed 
that the TLS-TLS coupling in glasses is dominated at 
low temperatures by an up-down transition (cf. Fig. 1). 
Such an interaction contains coherence because the up 
transition of one TLS inevitably has to be followed by 
a down transition of the coupled TLS. After rotating to 
the eigenbasis |0i), |li) of (n[) by 



iejE,)Sl-iX/E.,)Sl, 



(21) 
(22) 



where SI — \Qi){Oi\ — |li)(li|, such an interaction is gen- 
erated by the 



^'^3 T qi qj 



(23) 



part of (|T5|). The eigenstates |1) = \0i,0j), |2) = |0,, 1^), 
|3) = \li,Oj), |4) = |li,lj) of the Hamiltonian i?o,ij — 
(cf. Fig. 2), become mixed due 



In the up-down subspace 
|1,) and|l„0,) = |l,)®|0,) 



to the interaction term 
spanned by |0i, Ij) — \0i 
this pair coupling can effectively be described by a TLS- 
Hamiltonian of the type (|l|) with ct^ = |0i, lj)(Oi, lj| — 
|l„0,)(l„Oj|, CT, = |0„l,)(l„Oj| + |l„0,)(0„lj|, and 
pair asymmetry energy, pair tunneling frequency, and 
pair level splitting 



Cp = Ei — Ej, 
Ap = JijAiAj/2EiEj, 



(24) 



the single TLS in the glass. Burin and Kagant^l derived 
the following distribution function for the parameters of 
coherently coupled pairs 



P(2)(ep, Ap) = ^(p„fc^T)(PoC/o)^e(Ap 



A, 



■p.mm) ; 

(27) 



where Q{x), the unit step function, is included to empha- 
size that the distribution has a cutoff at small Ap. Note 
that the density of states of the pairs is linearly tem- 
perature dependent. Coherent coupling between pairs is 
destroyed if there is spontaneous decay during the up- 
down transition. For primary TLS with E w 2kBT this 
occurs with a rate (Q), which guided Burin and Kagan 
to estimate the lower bound (for secondary TLS formed 
from symmetric primary TLS) 



Ap,i„in = TT-aT 



Comparing (|2^) with (0) yields 

A{T) - (^Vl2) {P^Uo) keT. 



(28) 



(29) 



This looks encouraging; however, after a closer exam- 
inaticp. there are several inconsistencies. First, Maier 
et al.Ej could fit their data with a temperature inde- 
pendent parameter A. Burin and Kagan's theacy gives 
A(xT. The 10% variation, which Maier et aLB found 
between the 1 and 0.5 K data, is too weak to account for 
the linear temperature dependence of the theory. Sec- 
ond, putting in numbers, we find for PMMA at 1 K: 
A{1 K)/kB w 10""* K and Ap,min/fcB ~ 400 mK, which 



The eigenstates for the coherently coupled pair then read 



is inconsistent with the experimental values (13), (|lj) by 
several orders of magnitude. Indeed, at 1 K, TLS with 
Amin/fcs ^ 0.4 K can never be responsible for spectral 
diffusion on the time scale between 10'^ and 10^ s. How- 
ever, it is possible that very asymmetric primary TLS 
are responsible for long time spectral diffusion. For very 
asymmetric TLS, the estimate of Burin and Kagan (28) 
would be significantly reduced. Hence, though Burin and 
Kagan's theory predicts the measured time-dependence 
quite accurately, there arise severe inconsistencies in or- 
ders of magnitude and the temperature dependence upon 
applying their theory at 1 K. It should be mentioned that 
Burin and Kagan considered only the millikelvin regime, 
which avoids all these problems. 



1 + ep/Ep 



Jp/Ep 
2 ' 



0,), (25) 



<^p/Ep 



1 — Ep/ Ep 



|0„1,). (26) 



Clearly, the more asymmetric the pair is, i.e., the larger 
the energy offset Ei — Ej is, the more localized at one 
TLS is the pair excitation. 

Based on the distribution (||) of the single TLS tun- 
neling parameters and a uniform spatial distribution of 



III. THE KASSNER-SILBEY APPROACH FOR 
PRIMARY AND SECONDARY TLS 



Though the picture developed by Burin and Kagan is 
very appealing, it explains oply-jqualitatively the hole- 
burning data of Maier et al.EJ^t3. The question arises 
whether a microscopic calculation can yield testable pre- 
dictions based on the interacting TLS scheme of Burin 
and Kagan. We carry out calculations based on strong 



3 



coupling of TLS to phonons with deformation potential 
7 « 1 eV. This allows for a shifting of the rate distribu- 
tion towards longer times and indeed brings the "calcu- 
lated" value of A (see (29)) into closer agreement with ex- 
periment. Furthermore, the inclusion of strong coupling 
effects alone gives a reasonable fit to the experiment of 
Maier et al. for intermediate times (cf. Fig. 3(a)). 

Let us start with the usual spin-boson Hamiltonian for 
an ensemble of TLS interacting with phonons via a strain 
field. In the TLS-eigenbasis the Hamiltonian reads 



H = - 



with 



and 



i-yq 



f 1 



\ 2gVu!q 



1/2 



(30) 



(31) 



(32) 



According to this, the TLS become dressed with clouds 
of virtual phonons. As a result the coupling between the 
TLS will also be changed. The assumption tradition- 
ally made is that the dressed entities can be considered 
as weakly interacting. Hence, first order perturbation 
theory in the dressed states might be sufficient at low 
temperatures. i— . 

Based on this picture, Kassner and SilbeyllZl derived 
a iS'S^-interaction between TLS and a reduction of re- 
laxation rates for asymmetric TLS. Compared to Eq. 
(^ , they found well below the Debye-temperature of the 
glass, T < Gd, 



R 



where [cf. Eq. 



-G{e/E)^ n, 



G = (h/Snki 



D- 



(33) 



(34) 



It is a peculiar feature of their approach that symmetric 
TLS have no reduced rates, i.e., i?inax is still given by Eq. 
(^, and have zero interaction. It is this very fact which 
significantly changes the distribution function P^^^ {E, r) 
of single TLS energies E and dimensionless relaxation 
rates r = i?/i?max compared to the STM-result (cf. first 
term in Eq. (||)). According to (^3|) 



(35) 



(compared to the STM-result r ^ {A/E)'^). With this, 
Kassner and Silbey derived the new distribution 



P'-^\E,r) = 



Po 



2ru{r){l + G[l ~ u^{r)]}' 
where u{r) = e/E is the inverse function of 

r{u) = (I-u2)e"^"'. 



(36) 



(37) 



The result is a stretching of the distribution in r towards 
smaller values, i.e., P^^^{r) has an extended tail for such 
rates. The flaws and merits of the Kassner-Silbey way of 
handling strong-TLS-phonon-coupling effects have been 
discussed in Ref. |T^. An important point is that a SJ^Sl- 
interaction between TLS, as used in Burin and Kagan's 
approach, cannot be derived from this approach. 

In Appendix A we have generalized Kassner and Sil- 
bey's approach to include coherent coupling between 
pairs. Our procedure has been as follows. First, we 
have eliminated the diagonal {SV) and off-diagonal (5^) 
coupling in the Hamil toni an (pO ) by two separate uni- 
tary transformations (A.l) and (A. 6). Second, instead 
of continuing with the full transformed Hamiltonian, we 
have projected out the one-phonon coupling part in the 
dressed state basis. This has been achieved by expanding 
the generated phonon shift operators around their mean 
value up to the first leading term in c^. This generates 
one-phonon transition matrix elements in the four-level 
system of the pair and, more importantly, Debye- Waller 
factors which renormalize Burin and Kagan's coherent 
coupling term (23) as 



J,,A,e-<=(^-/^")'/2A,e-«(^^/^^)'/V(4i?,i?,)^^^^. 



Here, 



huJn 



(38) 



(39) 



which is equivalent to (|l^)-(|2^). The pair asymmetry 
energy and tunneling amplitude read 

ep = i?,e-^'«"?/2 - i?,e->^^G«?/2^ (40) 

Ap = JyA,e-'="'/2A,e-^"?/V2S^Sj-, (41) 

with W, = e"'^"'/2. One should note that for G > I the 
Debye- Waller factor g"^-";/^ is practically unity except 
for nearly symmetric TLS. Since we will be mainly in- 
terested in strongly asymmetric TLS, we will always use 
this simplification hereafter. In Appendix B we calculate 
the pair parameter distribution function. The result is: 



p(2)(ep,Ap)=p(2)(Ap)p(2)(ep), 



(42) 



with 



p(2)(Ap) = ^ e-'^/2/„2(G/4) e(Ap - Ap,n,in), (43) 



A2 

p 



4 



where Io{z) is a modified Bessel function, and 
with 



max ) 



(44) 



(45) 



1 + e~^"p 



log 



1 + e-P^- 



1 -|_ ^-P[<ip+E„ 



The parameters £"„ 



are cutoffs which are set by 



the requirement of stabihty for the coupled pairs. i?max 
is a cutoff in the energy of the primary TLS, which, as we 
will see below, is generally a function of the temperature. 
If i?max ^ ^bT, P^'^\ep) is the sum of a constant and a 
bell shaped curve centered around ep = with a width 
of O^ksT). It can be approximated by 

P(2) (ep) « A{T) [log 2 + (1 - log 2) eM-f3'el/9)] (46) 

for all practical purposes. If i?max ksT, the distribu- 
tion P^'^\ep) becomes independent of e^, 



P(2)(ep) ~ A{T)E^,^{T)/4kBT. 



(47) 



In the limit -Emax ^ ksT and GjCp/ksT <C 1, we find 
Burin and Kagan's result (p7|). This confirms that their 
model is valid at ultralow temperatures for nearly sym- 
metric TLS with weak coupling to phonons. In the limit 
G 3> 1, the asymptotic expansion of the modified Bessel 
function, Io{z) fa e^/\/2nz, [z 3> 1), provides a renor- 
malization of Pq by a factor 2/{ttG). 



IV. STABILITY ANALYSIS 

To determine -Emax and Apminj we first neglect the 
level broadening effect of the phonons. Then, -Bmax 
scales with the glass transition temperature, and Ap^min 
is set by the concentration Np of pairs in the glass. 
To determine Np^ we start with estimating the proba- 
bility for 2 TLS separated by a distance r to form a 
pair to P(r) « PoC^o/(87r/3)|rp. Then the concentra- 
tion of pairs in a shell \V — V /2,V + V /2] is given by 
Np Ki j^l"^ nP{r)dr, where n — P^ksT is the number of 
thermal TLS. This yields 



7Vp = (l/2)(PoT)(Po(7o)log3, 



(48) 



which in three dimensions is independent of the volume 
V due to the |r|~^-law of interaction. The maximum vol- 
ume K = (47r/3)|r*|'^ each of these pairs can occupy is 
given by 1/iVp which determines the minimum interac- 
tion energy Jmin = C/o/|r»p. This provides according to 
Eq. (H) 



1 TT 2 

Ap_min — — -^{PaUo) ksT 



(49) 



as a reasonable estimation of the lower bound. Indeed, 
for PMMA, Eq. (||) provides Ap^^ijke = 10"^ K 
which, at 1 K, corresponds to a maximum relaxation time 

Tmax = 10*^ to 10^ S. 

Let us now include decoherence effects by phonons. 
We will not provide a full discussion of the relaxation 
dynamics in the four-level system of the pair, but in- 
stead try to arg ue m ore p hysic ally. After two unitary 
trans forma tions (AT) and (A^), one finds the Hamilto- 
nian ( A.ll| ) as is pointed out in Appendix A. Here, only 
the term oc Sx'^ Sy allows relaxation within the up-down 
subspace {|0, 1), |1, 0)}. The relaxation mechanism is a 
fiip-flop process linked with the emission or absorption of 
a phonon. The relaxation rate scales with (Ap/Ep)'^: 



where 



d(p) _ „ d(p) 

l\, — 'pJl-maxi 



Ap 

Ep 



(50) 

(51) 
(52) 



and Xp — Ep/2kBT. We used that ap ~ a, as discussed 
in Appendix A. Note that there is no Debye- Waller fac- 
tor of the Kassner-Silbey type because of the coupling 
of the phonons to instead of to erf. The coherent 
coupling is destroyed by spontaneous decay during the 
up-down transition of primary TLS constituting a pair: 
for instance |0, 1),|1,0) |0,0). Based on this argu- 
ment, Burin and Kagan estimated Ap.min = haT^, where 
aT^ is the decay rate of symmetric, primary TLS with 
E = 2/cbT. Clearly, on a the time scale explored in 
the experiment of Maier et al., strongly asymmetric TLS 
dominate the hole width instead of symmetric ones. Note 
that in constructing P(^)(ep,Ap) (Eq. (42)), we inte- 
grated over values of the energy splittings of the primary 
TLS less than -E'rnax(T'), which have the correct initial 
population factors to insure the creation of stable pair 
excitations. We now investigate which primary TLS are 
able to guarantee stability of secondary TLS at 1 K, i.e., 
we ask whether the limit i?max(P) > ksT or < ksT ap- 
plies. We first study if secondary TLS may be formed 
from primary TLS with thermal splitting, E ^ O^ksT), 
at 1 K. We need to satisfy three conditions. First, we 
require 



Ep/h > R{E = 2kBT), 



(53) 



where R{E = 2kBT) = raT^<l/t is the relaxation rate 
of a primary thermal TLS which has not yet decayed at 
t, i.e., for which (cf. Eq. (||)) 



1/r > aT^t 



10i3 - iQie 



(54) 



with t — 10^ — 10® s. These two requirements guar- 
antee that the secondary TLS is coherent on the time 



5 



scale where deviations from logarithmic spectral diffu- 
sion is seen. Furthermore, if secondary TLS are respon- 
sible for the spectral line broadening after 10^ s, the pair 
rate must satisfy R^pH — 1. This provides the relation 
(Ap/Ep)'^ — 1 / aT^tXp coth x„ . Using this relation, mul- 
tiplying Eq. ( p3[) by ^/Ep and squaring it, one can now 
easily check that Eq. ( ^ ) always implies Eq. (|53|). 

With respect to the stability of secondary TLS at 1 K, 
we conclude: First, from the condition (54) we deduce 
the criterion r <C 1, i.e, the primary TLS must be very 
asymmetric. Second, upon replacing Ap = Jr [cf. Eq. 

(p])] and noting that Xp/ \J x^ coth Xp w 1 for < < 1, 

we find the criterion 



J 



2kBT 



> VaTH ^ 10^ 



10« 



(55) 



Hence, according to J = Uo/d^ and Uq w IQ-* K A^, 

the relevant TLS are separated by a distance d ^ 1 A. 
Note that these estimates are highly approximate. In 
particular, we have approximated the rate at which co- 
herence in a secondary TLS is destroyed by the relaxation 
rate of a primary TLS, without consideration of the true 
rates that govern the coupled four level system. How- 
ever, we may still conclude that at 1 K, it is unlikely 
that primary TLS with thermal splitting, E ~ OlksT), 
can form pairs, unless they sit extremely close to each 
other. This fact suggests a natural cutoff, -Emax(r), in 
the energy splittings of the primary TLS that comprise 
the secondary TLS existing at 1 K. This cutoff is moti- 
vated by the above stability criteria, and can change the 
temperature dependence of P^^^ (e^) in Eq. ( |4^ ) and of 
the hole width in general. We will return to this point 
in the next section. According to this, the distribution 
P'-^\Ep,rp) reads 



P^'HEp,rp) = P<^'\rp)P<^'HEp,fp) 



(56) 



where 



Poe-J/g(G/4) 
2rp 1 Tp 



0(rp - rp^inin), (57) 



and 



P^'HEp,rp) 



A{T) 



(58) 



X lo, 



log 



l + e 



-0E„ 



(Ep ^ -Ep), (59) 



with < Ep < i^max(T), rp = ^1 - r, and A{T) given 
by Eq. M). 



V. COMPARISON WITH EXPERIMENT 

We now calculate the broadening of a spectral hole 
in the inhomogeneous line due to spectral diffusion in- 
duced by sijQgkj TLS and pairs. According to the stan- 
dard theoryLTtiJ, the experimentally observed line broad- 
ening, Ar(t) = r{t) — r{to), can be written as 



Ar(i) = ^(C) 



2 E 
dE sech — — — x 
2fcsT 



[u{r)P^^^ {E, r) + 



X I dr 



(60) 



where P^'^\E,r) Sind P^'^HE,r) is given by (H) and 
respectively. To calculate these integrals, we replace the 
last factor by the step function, which restricts the r in- 
tegration to the interval [l/i?max^j l/-Rmax^o]- This gives 



Ar(t) = Ar(i)(i)-t-Ar(2)(t), 



where 



ArW(0 = J^(C)PofcsTlog 



3h 



l-u'^{l/aTHo) 
l-u^{l/aTH) 



(61) 



(62) 



is the contribution of the primary TLS to the spectral 
hole width. To calculate the pair contribution Ar^^\t), 
we consider the limits E^i^^{T) 3> kgT and <^ k^T, sep- 
arately. 



A. The limit E^^^{T) > ksT 

Remember, in this limit coherently coupled pairs can 
be built by thermal primary TLS. We find from Eq. (|6C" 

Ar(^H^) = ^e-^/%^(G/4)(C)Po X 



A{T) (Vt - 



(63) 



Putting Eqs. (pll)-(|6^) together provides an equation 
which is very similar to the result (pi]). However, this 
equation is based on the microscopic picture of phonon- 
mediated TLS-TLS interaction. There are two regimes 
where the function u{l/aT^t) can be determined an- 
alytically: (i) short-time limit: aT^t e*^ which 
gives u'^{l/aT^t) w YfGlog{aT^t); (ii) longtime limit: 
aTH > which gives u'^{l/aTH) « 1 - e^/oTH. In 
Ref. 12(b), Maier and Haarer have fit their data with 
only the first term in (|6l|). By numerical inversion of 
Eq. (^rf), they could find good agreement on intermedi- 
ate time scales up to 200 — 300 min. We have illustrated 
this in Fig. 3(a). The values for the fit parameters G 
and Po(C) are 32 and 6 x 10~^, respectively. With this 
value for G, the crossover from the short- to the longtime 
behavior happens after ca. 150 min; Eqs. (|l|) und (11) 
yield a Debye-temperature of 9d = 108 K, which is in 
reasonable agreement with literature data for PMMA. It 
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is in part due to the success of the Kassner-Silbey theory 
at intermediate times that we have adopted the strong 
couphng approach as our dynamical starting point. For 
times larger than 300 min, Maier and Haarer attributed 
the deviations of the theory from the experimental data 
the contribution of interacting TLS. If this is true, we 
should find agreement between theory and experiment 
when including the second term i n (|6l| ). In Fig. 3(b) we 
have plotted AT{t), Eqs. (p|)-(|63|), together with the 
experimental data of Ref. |12| for 1 K (upper curve) and 
0.5 K (lower curve). We used the same value for G as 
in the previous plot, and have optimized the TLS-TLS 
coupling parameter, PqUq, and the TLS-chromophore 
coupling parameter, Pq{C), to find best agreement for 
the 0.5 K data. The result is PqUq = 2.5 x lO"*' and 
Po{C) = 4 X 10~^. The upper curve shows the predic- 
tion of the theory for the 1 K data. The parameter value 
for PqUq, which fits the 0.5 K data^jJs by a factor 25- 
35 smaller than the literature value.H However, we see 
that a superposition of a logt- and a t^-^-term can be 
interpreted as an effective t^ '^^-power law on the exper- 
imental time scale as seen by Maier et al. Hence, under 
the assumption that primary TLS with thermal energy 
splitting can form stable secondary TLS at 1 K, we find 
AF^^^ cx T^/^, giving a temperature dependence that is 
too strong compared with the experimental observation 
which indicates AF*^^^ cx T'^/^, at least in fitting the data 
at the two temperatures, 0.5 K and 1 K. This can clearly 
be seen by the prediction of the theory for the data in 
Fig. 3(b). 



B. The limit _Bn,ax(T) <^ kgT 

We have noted in the previous section that it is very 
unlikely that thermal primary TLS can form stable pairs 
at 1 K. If we impose a minimum separation distance be- 
tween primary TLS (say 5 A) , then a natural energy cut- 
off enters due to the stability requirements outlined in 
the previous section. If we assume that this cutoff -Emax 
satisfies the condition 



and 



£^max(T) « fcsT, 



(64) 



then we can estimate £^max(T') from the conditions con- 
sidered in Sec. IV. Specifically, from Rt ^ 1, R^pH ~ 1, 

Xp/ ^ x-^ coth Xp ^ 1, and r = Ap/ J, 



2kBT 



coth 



2kBT 



J IT 



2kBTVaTH' 



(65) 



where J^ax = C/o/rfmin- This shows that £^max(r) 
T~^/^, which significantly alters the temperature depen- 
dence of the hole width. If the above condition is met, 
then 



AF(2)(t) ^ ^e-«/2/o(G/4) {C)Po {Vt - VUi) x 



(67) 



P^^HEp) « A{T)E^a^iT)/4kBT ^ T-^l\ 



(66) 



In the regime where Emax{T) <C fcsT, this expression is 
essentially temperature independent. Thus, the stability 
requirements imply an interesting thermal breakup of the 
secondary TLS. 



VI. DISCUSSION 

For very low temperatures (i?max(^^) ^ ^bT), the ma- 
jority of asymmetric TLS are stable even if constructed 
from asymmetric primary TLS that have energy split- 
tings on the order of fcsT, leading to the strong tem- 
perature dependence depicted in Fig. 3(b). We have 
argued that at higher temperatures, a crossover should 
occur where the temperature dependence should become 
weaker, as secondary TLS become less stable. Therefore, 
it is qualitatively consistent with this picture that the ob- 
served temperature dependence at ~ 1 K is weaker than 
that shown in Fig. 3(b). Such arguments require further 
"slowing down" of the temperature|-d|ependence. This is 
indeed seen experimentally at 2 K.l3 Hannig et al. are 
able to fit their data with essentially the same value of 
A at 0.5 K and 1 K, and a value of 0.3 x A{T = 0.5 K) 
at 2 K. Eventually, the entire TLS picture should break 
down at some temperature in the range 1-10 K. A clear 
test of the validity of this picture would be to see if the 
stronger temperature dependence emerges at lower tem- 
peratures. Such experiments are difficult to conduct at 
ultralow temperatures, due to slow equilibration effects. 

We may also note that the discrepancy between the 
"derived" and literature values of PqUq can be explained 
by the thermal breakup of secondary TLS demanded by 
stability criteria. As the temperature is raised, the pri- 
mary TLS must lie closer together to form a stable sec- 
ondary TLS at long times. The probability of finding 
two neighboring primary TLS a very close distance apart 
is small, which effectively decreases the "derived" value 
of PqUq- Note that there is also a reduction in AF(^)(t) 
due to the factor {E^a^/2kBT f (see Eq. (|6^)). Thus, 
at least qualitatively, it is indeed possible that a picture 
based on coupled TLS can account for all the properties 
currently observed in longtime spectral diffusion experi- 
ments. 

So far we have shown the pros and cons of applying 
the idea of coupled pairs to hole-burning experiments 
on ultralong time scales at 1 K. Here we propose an 
interesting experiment which, although difficult to per- 
form, would provide a conclusive test of the coupled TLS 
hypothesis. A crucial observation is that the exponent 
of the power law depends on the dimensionality of the 
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glassy probe. For example, the exponent would be sys- 
tematically smaller if the primary TLS were confined 
to (quasi)-two dimensions, while the interaction between 
them still varies as l/|rp- This would result in a distri- 
bution function -P(Ap) oc Ap^"^ (from P(Ap) c>c P{J) — 
P(|r|)|r|(|d|r|/dJ|) and J cx |r|-3^ P(|r|) = constant) 
and therefore in a t^/'^-power law. The experiment we 
suggest has already been performed, albeit not not for 
the purpose that we discuss here and not on a time scale 
upf-tjo 10^ s. The hole-burning experiment by Orrit et 
al.u on an ionic dye in a Langmuir-Blodgett monolayer 
lying on a three-dimensional substrate is an experiment 
of the type we proposed above — the TLS dynamics is 
restricted to two dimensions whereas sound waves and 
strain are not affected by the interface between the-amor- 
phous layer and the bulk. Indeed, Orrit et alE^ have 
observed spectraLdiffusion which could be explained by 
Pack and Fayero using the standard tunneling model. 
From this perspective it seems promising to extend Or- 
rit 's experiment to longer times and look whether a power 
law weaker than the three-dimensional result t'^-^ could 
be observed. 

There are, however, reasons to be skeptical of the pic- 
ture we have outlined. A number of approximations have 
been invoked that render only a semiquantitative descrip- 
tion of the experimental results. These approximations 
include the reduction of the primary TLS to an effective 
secondary TLS, and the use of relaxation rates for uncou- 
pled TLS to discuss decoherence effects for coupled sec- 
ondary TLS. These approximations, especially the first 
one, may not be valid at temperatures near 1 K. Clearly, 
one has to think about alternative explanations. A sim- 
ple idea would be to attribute the deviation from the 
standard log i-behavior to a nonequilibrium state of those 
TLS which relaxes on these long tinM-|Scales. Indeed, re- 
cent experiments by Friedrich et al.E3 have proven that 
spectral diff'usion in glasses under nonequilibrium condi- 
tions results in a nonlogarithmic time evolution of the 
hole width. However, the data in Maier's experiment 
were obtained after letting the sample relax at the mea- 
surement temperature for a longer period than the later 
data recording period. Hence, one expects that all re- 
laxation processes shorter than this waiting time occur 
under equilibrium conditions. 

It is interesting to compare the glass results with 
equivalent experiments on proteins. Hole-burning exper- 
iments on proteins show almost no hole broadening up 
to three houra^oUowed by a nonlogarithmic hole broad- 
ening regimeO A new type of temperature-cycle hole- 
burning experimenta2a leads to the conclusion that the 
excess broadening of the hole in the protein cannot be 
interpreted in the framework of the STM. Hence, one 
might speculate that in both glasses and proteins, the in- 
teraction of TLS becomes important at long time scales. 
An alternative conclusion avoids the notion of interact- 
ing TLS altogether. Perhaps the energy landscape of 
glasses is not too diflFerent from that of proteins, and also 



shows features of organization at high barriers. Recall 
that the experiments of Maier et al. were carried out in 
PMMA, which is a polymer glass. Polymer glasses may 
be expected to have conformations similar to proteins. 
For example, such "conformational" dynamics may in- 
volve side chain motions. The physical picture is that 
the energy landscape comprises high barriers in addition 
to constantly distributed lower barriers within each of its 
basins. The algebraic behavior then results from tunnel- 
ing through those high barriers, which have to be dis- 
tributed around a "typical" value Vq. This value has 
to be sufficiently high in order that the onset of the al- 
gebraic behavior occurs only after 3 h. We give some 
details on these ideas in Appendix C; more can be found 
in Ref. ^ The model predicts a temperature and time- 
dependence of the hole width with an exponent which 
is slightly weaker than T^/^ and t^^^, respectively, and a 
slowly decreasing function of T and t. We note that if in- 
deed specific polymer dynamics (like side chain motion) 
are responsible for deviations from logarithmic spectral 
diffusion, perhaps a deuterated sample may show differ- 
ent hole broadening behavior. 

One is tempted to speculate that these "non-STM like" 
high barriers arisjg-from the presence of the chromophore 
in the glassy hostEII, because they have not been-observed 
in sound attenuation experiments up to 100 K.E£lt3 Note 
however, that these experiments were not performed on 
polymer glasses, and that sound attenuation experiments 
on PMMAr-a|nd PS indeed show an additional peak above 
ca. 50 K.e£| Interestingly, doping a network glass with 
OH-impurities leads to the same observation^^, which 
conceivably supports the importants of side chain mo- 
tion in polymer glasses. Furthermore, in contrast to the 
glass, the protein probe is not doped by a dye. Instead 
part of the protein is chemically changed in order to serve 
as a chromophore. Thus, it seems that such deviations 
from standard spectral diffusion behavior are not due to 
the inclusion of the chromophore into the sample. 

One should also note nonlogarithmic line broadening 
is the typical case because a logarithmic time dependence 
occurs only in case of a 1/ A"-distribution for the singu- 
lar case n = 1. Without invoking the physical reasons 
for deviations from the standard 1/A distribution func- 
tion introduced by Anderson, Halperin and Varma, and 
of Phillips, we may say that for small values of A (corre- 
sponding to long times), the distribution of barriers is not 
really flat, but instead a smoothly varying function of the 
parameter A (see Appendix C). That the distribution of 
barriers in such a model turns out to be a log-normal dis- 
tribution shows a striking similarity to general systems 
exhibiting 1//" noiscEd 



VII. CONCLUSIONS 

In this paper we have analyzed the consistency of the 
conjecture that coupled pairs of TLS dominate spectral 
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diffusion on ultralong time scales. Because the pairs 
distribution P*^^-'(ep, Ap) is correlated with the distribu- 
tion of the primary "STM-like" TLS, we have in a sense 
pushed the STM as far as possible by looking at these 
new low energy excitations. We believe that this is an 
important step, which has to be done before trying to find 
another extension of the STM for every new experiment. 
The questions was whether they also apply to the 1 K 
regime, and, in particular, whether they can provide an 
explanation of the longtime hole-burning experiment of 
Maier et al. We find that a picture based on the idea of 
interacting tunneling systems seems consistent with the 
experimental data, although we are unable to fit certain 
aspects of the experiment, such as the temperature de- 
pendence, quantitatively. Also, alternative explanations 
have been presented. These models are, at the moment, 
at least as speculative as the scenario of coherently cou- 
pled pairs. For this reason, we have discussed some pos- 
sible experimental tests of the theoretical models we have 
presented. More theory and experiments have to be done 
to finally evaluate the role of coupled TLS in glasses, and 
to understand the origin of nonlogarithmic hole broaden- 
ing in glasses. 
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APPENDIX A: 

In this Appendix we generalize Kassner and Silbey's 
approach to include coherent coupling between pairs. For 
this goal we first apply the unitary transformation 



where 



Ui ~ exp 



]i2u,cyhL,,)ib,-bl„)Sl/2 



l,q 



(A.l) 



to the Hamiltonian ( ^ ) which eliminates the diagonal 
coupling of phonons to SI- The transformed Hamilta- 
nian H U^^HUi reads with Si = (1/2)(S'J ± iS^H 

H^-^J2^^^l + T.^''i^% (A.2) 

i q 

- J2 + blq) [BISX + B\S^_) 
hi 

- jY.-'u {mjSlS;-2uiUj(B'_Sl + BXSl)Sl), 



q -q 



huJn 



(A.3) 



and 



exp 



±Y.{2u,djnu;,){b,^b\) 



(A.4) 



The S!^ are the usual phonon shift operators. Instead 
of proceeding with the full transformed Hamiltonian, 
we project out the one-phonon fluctuations around the 
shifted harmonic oscillators coordinates. This is achieved 
by expanding = (S^) + (S^ - (B^)) = (S^) ± 
where 



-G(e./iJ.)V2 



(A.5) 



and neglecting two-phonon terms. Applying Fermi's 
Golden Rule to the remaining one-phonon term directly 
yields the rate (|^. A coherent coupling between pairs 
is generated by eliminating this term by a second unitary 
transformation 



U2 = exp 



^(2u,(Si)c;Ac.,) {b,-bl^)Sl/2 



l,q 
PiSi/2 



(A.6) 



This yields: 

i q 

-^E-^^ {uiU,SlSi+u,{BV)u,{Bi)SlSi 

id 

-2u,u,(^{Bl)Sl-mSl)si}, (A.7) 

where 



SI = SI cosh (f)j — iSy sinh (f)j , 
Sy = Sy cosh (pj + iSl sinh (pj . 



(A.8) 
(A.9) 



Expanding again around the one-phonon fluctuations 
around the shifted harmonic oscillator coordinate by re- 
placing Si « SiiDi) ~ 10,5^ and 5^ « + i0,5|, 

where 



(A.IO) 



finds 
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H 



■i q i 

-J2j., [u,{D^^)u,{D^^)SlSi 



2iu,(Dl)uADi)ipjSlSi 



(A.ll) 

The only operators which act nontriviaUy in the up-down 
subspace{|0, l),|l,0)}areS'2(8)l, 105'^, Sx<»Sx, Sx<»Sy, 
and Sy 1^ Sx- If we define new pseudo-spin operators for 
pairs 



aP = |0„l,>(l^,0,| + |l„0,)(0„l 



aP = -i|0„lj)(l^,0,|-f i|l„0,)(0„l,|, 
aP = |0„l,)(0„l,|-|l»,Oj)(l.,0,|, 



(A.12) 
(A.13) 
(A.14) 



we can project the quoted operators onto the up-down 
subspace, which yields Sz 1 <yP, 1 ® Sz — cr^, 



0l,Sx^ 



1 Sy 



-qP and Sy ( 



project the Hamiltonian (A.ll) accordingly, we find for 
each pair 



aP. If we 



pair 



A.j-|.f + .A,.;i:|4(i,-.l,) 

g ^ 

^J2hu;,b%, (A.15) 



with 



Jij 



A,;A 



2EiEj 



EAD\ 



cP^c\ {ejE,)-cl{e,/E,). 



(A.16) 

(A.17) 
(A.18) 



Apply ing Fermi's Golden Rule to the interaction term in 
( A.lSj ) provides a one-phonon relaxation rate 



Rp = aAlEp/{2kBf coth{Ep/2kBT). 



(A.19) 



Here, we have used that the oscillating term in c^c^* 
is small so that pairs and single TLS relax approxi- 
mately with the same coupling parameter a provided 
that ei/Ei « 1. Note that because of the coupling of 
phonons to (instead of a^), there is no diagonal cou- 
pling term cx as in ( |30| ) after rotation to the pair-TLS 
eigenbasis. 



APPENDIX B: 

In this Appendix we calculate the distribution func- 
tion for the pair parameters Ap and Cp of Eqs. ( |40| ) 



and (^l|) . We note that the single-TLS distribution func- 
tion for u = A/E and E = A? + reads P{E,u) = 
Po/{uV^ — u^)- The distribution function for the pair 
parameter is then given by 



P(2)fe A) - [ f 



dEo 



I dui [ du2 [ dJP{Ei,Ui)P{E2,U2)P{J) 

Jo Jo J 



X - {Siep -Ei+ E2) + S{€p + Ei- E2)} 

■G{\-u\)l1-G{\~u\)l2 



X (5 ^Ap — \ Ju\u2e 



(B.l) 



The factors (l-l-e^'''^)"^ account for the thermal occupa- 
tion of the primary TLS. If the TLS are homogeneously 
distributed in the glass and interact via a dipolar cou- 
pling, J = J7o/|rp, the distribution function of J reads 
for a three-dimensional probe 



P{J) 



47rJ7o 



Using now that 



du^-^^ = e-G/%(G/4), 



(B.2) 



(B.3) 



where /o(^) is a modified Bessel function, and 



dx 



1 



log 



1 + y 



log 



1 + M 

1 + My, 



with X — e , y — e ^^""^^^ and /i = 
finds Eqs. (||)-(|5l). 



one easily 



APPENDIX C: 

In this Appendix we show how a modified tunneling 
model which comprises an additional Gaussian distribu- 
tion in the tunneling parameter A centered around the 
mean value Aq ^ 1 can provide an algebraic line broad- 
ening. The new distribution reads with A = hujo e^^: 



P(e, A) = Po + 



PqA 



(C.l) 



where A is a dimensionless constant, and A > \m\i\{E), 
Amin(-E) = \og[hujQ / E) . A similar niodel has been used 
previously|_by Jankowiak and SmallEZI and by Ziindars 
and Fayero to discuss 3-pulse photon echo data.E3'Lj A 
combination of both terms is needed in order that the 
onset of the algebraic line broadening does not occur too 
early. In the following it will turn out sufficient to reduce 
the number of parameters by setting 



= Ao. 



(C.2) 
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The ensuing distribution in A then reads 

1 S(?it^Q)l-(l/2Ao)log(Ra;o/A) 



/^2-(l/2Ao)log(;i(^o/A) 



, (C.3) 



where B = /^27rAo. This distribution has to 

be compared with (R). Defining the maximum relax- 
ation rate at given E hy R = e-2(^-^™(^))i?„iax(-E), 
the distribution function in relaxation rates, R, and TLS- 
energies, E, reads 



{e/E)PiE,R)^^ 



1 

R 



B R''^^'^ 



where 



i^iR) 



8Xq 



log(i?max/i?). 



(C.4) 



(C.5) 



From this expression, it is obvious that an algebraic line 
broadening Ar(i) cx with an exponent /i < 0.5 occurs 
for Ao ^ \og{Raiii^{T)t) . The exact calculation reveals 



3h 
A < erfc 



Xo~{l/2)log{KTt) \ 
V2X^ J 
Ao-(l/2)log(ifTto) 



(C.6) 



where erfc(a;) is the complementary error function, and 
K = {nujo/kefa. With htoo/kB = 0(1 K), one finds 



10^° K~^s- 



For Ao » log(-ftrrt) this gives an 



K 

algebraic line growth 

Ar(t) OC Sr(i^Tt)(l/2)-(l/8Ao)log(KTt)_ (C.7) 

More details and a comparison to experiment can be 
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FIGURE CAPTIONS 

FIG. 1: Coherent up-down coupling between two TLS. 

FIG. 2: Energy levels and eigenstates of the Hamiltonian 
Ho^.j = -{\/2){E,Sl + EjSi), {E, > Ej). Framed 
are the up-down states which build a basis for the 
secondary TLS. 

FIG. 3: (a) Hole broadening in PMMA at 1 K (upper 
curve) and 0.5 K (lower curve) compared with the 
Kassner-Silbey theory for single TLS. The exper- 
imental data are from Ref. 12, The two solid 



lines are Ar(i)(i), Eq. ^ for 1 (upper solid 
curve) and 0.5 K (lower solid curve) with G = 32 
(915 = 108 K) and Pq{C) =6x10-^. (b) Same fig- 
ure as above but now with the contribution of both 
the single and pair TLS [logarithmic and algebraic 
part of Eq. ^)]. The parameters Po(C) =4x10"^ 
and PqUo — 2.5 x 10~^ have been optimized to find 
best agreement with the 0.5 K data; G = 32 has 
been kept fixed. Comparison with the 1 K data 
clearly shows that the temperature dependence as 
predicted by Eqs. (29) and ( |6l| ) is too strong. 
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FIG. 2: Neu et al. 
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FIG. 3(a):Neuetal. 
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FIG. 3(b): Neu et al. 
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